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The exchange coupling and spin waves in metallic magnets: removal of the long-wave 
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Abstract 

A well-known connection between the magnetic susceptibility and the effective exchange parameter is analyzed. It is shown 
that all current computational schemes use a long-wave approximation which is suitable only for localized moments systems. A 
corresponding smallness parameter is derived in real and recipical space. General 'inverse susceptibility' approach is combined 
with a multiple scattering theory and applied for the studies of elementary metals Fe, Ni and Gd. It is shown that the proposed 
approach significantly improves the description of the exchange coupling between nearest atoms, the spin wave spectrum at 
large wave vectors for the itinerant degrees of freedom in metallic magnets. A consistent usage of this method leads to the 
proper description of magnetic short-range order effects. 
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The exchange coupling is one of the most important fundamental interactions in magnets. Knowledge of the 
parameters of this interaction facilitates the description of numerous properties of magnets and provides a critical test 
of the applicability of different models. In addition, the current applied research in the area of magneto-recording, 
spintronics and magnetotransport phenomena requires the knowledge of the effective exchange coupling between 
atoms or layers Therefore reliable calculations of the parameters of this coupling are especially important both 
from fundamental and applied points of view. Many computational expressions for magnets with different magnetic 
orderings were proposed^, 0, El IE IE IE However, most of them were derived only for localized spin systems 
at equilibrium. In this treatment the atomic magnetization is directly proportional to the effective 'exchange' field 
(localized approach) . In real non-equilibrium state these directions are always different due to coexistence of localized 
and itinerant degrees of freedom. Any model assumption about this non locality restricts certain degrees of freedom 
and affects the spin-spin correlation function and hence the dynamic and thermal properties. For instance, the 
local approach implemented in Ref.jE] directly leads to the concept of a strong short range magnetic order which 
by no means is a universal behavior of magnets. It is desirable that both short-range (large q) and long-range 
(small q) fluctuations be present in the general theory in a consistent fashion at any finite temperatures, whereas 
the perturbation theory 0, IE IE IE 13 is used only at small temperatures (spin wave (SW) stiffness calculations). In 
addition, for weak magnets or magnets with some weakly magnetized regions the non-spherical parts of the effective 
magnetic field must be included. 

In this paper starting from the general principles I provide a critical analysis of the approximations used in the theory 
of exchange coupling, and discuss an opportunity to calculate the parameters of such coupling rigorously without 
specific assumptions about the range of magnetic order or any approximations about the form of magnetization 
density. 

In a nonuniformly magnetized medium, the total energy depends on the relative orientation of vectors m(r) at 
different points. The stability of the magnetic state with respect to the variation of magnetization density is determined 
by the tensor of the second derivatives 

5^ E 

J af3 (r,r') = - Sma ^ Smp ^ \m=m (1) 

where a — x, y, z and mo is the magnetization of any stationary state. In general, unlike the parameters of Heisenberg 
model, the parameters J a p(r,r') depend on the magnetic configuration and are not necessarily related to the total 
energy differences between different magnetic phases. The quantity J a p is a model quantity and the definition Q is 
not unique. 

To calculate J a p (r, r') one has to know the variation over the spin density. The real physical perturbation which 
can induce such a variation is the external magnetic field. Thus, we have to employ the ideas of linear response and 
calculate the dynamic magnetic susceptibility 

(DMS)0. 

In this case the variation of magnetization is written as 



5m = x^ext, (2) 

where x is the (non-local) DMS which can be presented as 

^ dm _ S 2 E _ ~. 

-X f _ , s f _ X = XJ (r, r )x, (3) 



SB ext 5m(r)<5m(r')' 



Therefore J = %~ .This definition of the effective exchange is nothing but the well-known connection between the 
inverse static susceptibility and the parameter $ 2 in the Landau theory of phase transitions (T.Moriya also uses the 
inverse susceptibility in the Heisenberg model; see Ref.0])- The goal of this paper is to consider the relationship 
between this definition and those currently used in magnetism research. Also the corresponding modification of the 
adiabatic SW spectra will be analyzed. 

Let us generalize this concept and introduce the frequency dependent exchange coupling J(r,r",w) as 

dr"j(T,T",co)x(v",v',uj)=S(T-r'). (4) 

The parameter J(r, r',u) (below we will omit matrix notations where it is not important) is the full exchange 
containing all enhancement effects. Some approximate expressions for J(q,w) in a long wave approximation have 
been obtained in Ref.||. Unfortunately, the exact calculation of DMS in the case of interacting particles is not 
possible at the moment and any practical evaluations of J(v,r',u>) require certain simplifying assumptions. 
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Using linear response for interacting electrons for the inverse DMS one can decompose the total exchange as follows: 



J(r, r', u) = J°(r, r', u) - I xe (r, r', w) 



(5) 



Here J°(r, r', uS) is the 'bare' exchange coupling (obtained, for instance, from the Kohn-Sham wave functions) whereas 
I xc (r, r', a;) is the exchange and correlation tensor that connects B xc to the induced magnetization B xc (r,a;) = 
/ dr'I xc (r,r',uj)m(r',uj) . 

Eq[5] represents a very convenient way to add enhancement effects to the exchange coupling parameter. In the 
absence of an explicit form for I xc from many-body theory, Eq0 must be considered as a formal definition of the 
frequency dependent Stoner parameter I xc (y,y' ,uj). It becomes practical only when the translational invariance of a 
perfect crystal J(r,r',ui) — J (r — t',u) and the condition I xc {r, r', uS) = I xc (r — r',0) = I are assumed. Then the 
condition of the SW existence is written as 



J°(q,w) = 4 



(6) 



One has to emphasize that Eq. (J5J is suitable for full-potential band structure treatment and requires a complete basis 
set for its computation for the real materials. Also, below we will consider only the ferromagnetic case and will limit 
our consideration to analyzing the effects related to x H component. For arbitrary magnetic ordering, the matrix 
x(r, r', a;) is a 4x4 matrix and its x" 1 component can not be easily separated from \ zz an d charge components^. 

To illustrate how the parameter J enters the spin dynamics (SD) equation of motion (EOM), let us consider the 
adiabatic limit of EqJSJ In this case a following density functional torque equation is valid 



dm (r, t) 

dt 



: 7tn(r,t) x B(r) 



(7) 



where B the total field acting on the electronic spin at point r and 7 is gyromagnetic ratio. Not all components of m 
are independent in Eq[7| because of the relation m-dm/dt = 0. 

For the adiabatic SD Eq.Q can be solved for several simplest magnetic orderings by linearizing this EOM analo- 
gously to the rigid spin approximation (RSA) case 10] . For instance, for ferromagnet (FM) we obtain 



(8) 



where x q is the Fourier transform of the 'bare' static susceptibility %(r — r',u> = 0) in the case of a periodic system. 
Kq is a 'kinetic gradient' matrix which takes into account the difference in spatial gradients of wave functions for the 
different spin direction^]. This result for the adiabatic SW spectra is obtained here in the non-local case without 
using RSA for the shape of magnetic perturbation. For more complicated magnetic structures (spirals, and etc) the 
corresponding dispersion laws may be obtained using a rotating coordinate system. The above results were obtained 
assuming that the decay of SW is absent (adiabatic SD); otherwise these equations are exact (in linear regime) and 
treat both short and long wavelength scales on an equal footing. If a magnetic field of relativistic origin is included, 
then in Eq. (JHJl one can add the energy of the free precession (rotation of the magnetization of the whole crystal). 

To analyze the relationship between the dispersion law (JBJ) and thatjij commonly used in the DFT, ui^ — 
ml [xo — X q ] I> we assume that the ratio 



A = TrA = Tr ( Xq - Xo) X^ 1 » u'jml 



(9) 



is small (long- wave approximation) and Kq — 1 (effective bandwidths for the different spin directions are equal). 
Then, by expanding Eq|S]over the parameter A, we obtain the desired result: 



£ q = m(Xq 1 -Xo 1 )= m X Q 



w q [ 1 - 5q ( 7m 



x a(i-a) 

- mXo 1 (xo ~ X q ) Xo 



(10) 



m (4 - J q ) 



or 



Jq — Jo 



1 + A 









= Jo 



1 + A + A 2 



(11) 
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where the susceptibility has a matrix structure and 



4 = Jo (l + A) = Xo -x q Xo 1 = IX«I = -'^bS^ 7 (12) 

is the matrix of the exchange parameter in the local (long- wave) approximation. Ea. ll2|) is the widely accepted 
definition of the exchange coupling parameter which was introduced in Ref.0>0,Q and used in different modifications 
in Ref . 0, 0, 0, Q ■ EqlHland the relations (|f t)lf f \ arc the main results of this paper. Below we will apply this result to 
the calculations of the exchange parameters and adiabatic SW spectra in Fe, Ni, and Gd using the multiple scattering 
technique. 

From the computational point of view two alternatives can be used. First one is a calculation of the total energy of 
the system with an external field included and the second one is to use of so-called 'local force' theorem which does 
not require self-consistency and the total energy calculations. 

In Ref.0 a general 'local force' theorem was formulated for magnetic perturbations and a convenient multiple 
scattering technique was suggested for calculations of \ xc = —d 2 E/dH xc d'B xc , where a role of B IC belongs to one-site 
scattering matrix t . We stress that this theorem can be used for the calculations of both the magnetic susceptibility 
X H and the effective exchange J = —d 2 E/dmdm without any modifications of the original theorem. This statement 
can be also seen as a trivial result of the well known linear response theory. 

The usual assumption of weak enhancement (see EqEJ has always been used previously. This is a reasonable 
approximation for the SW spectrum, because if we assume usual I x Jsi) = !■> then the important equality J(q) — J(0) = 
J° (q) — J° (0) is valid. As a consequence, the spectrum of elementary excitations is not affected by exchange-correlation 
enhancement effects. So, the definition i|12|) is directly related to several very strong approximations: RSA, smallness 
of SW dispersion compared to the effective exchange splitting (atomic limit), and small dispersion of I xc (q,w). 

From the formal point of view, the RSA was directly contained in the way that the local force theorem was previously 
used. While the exact formulation of this theorem contains the first functional variation of the total energy, in RSA 
one has to use the corresponding first derivative 

SE 9E (13) 



5m (Rj + r) dra.i 

with all gradient terms (large q vectors) omitted. Here i is index of the atomic site. 

The second approximation (x q — Xo) Xo 1 << 1 also removes certain degrees of freedom with large q which are 
present in the Brillouin zone. Analogous one site approximations have been employed in modern dynamic mean field 
techniques^], self-consistent GW scheme^ and the coherent potential approximation^. Such a restriction affects 
the spin-spin correlation function and in turn the dynamic and thermal properties of magnets. It is true that (as one 
can easily obtain from Ea lll(l . such a long wave approximation can safely be used for SW stiffness calculations. 

The term I xc (q) is usually short-ranged, which coincides with the corresponding smallness criteria in real space 
Xij/Xu << 1 (where Xij is a Fourier transform of x q and Xu is the on-site susceptibility), hence the local approximation 
-fxe(q) = I niay affect the exchange between the nearest atoms only. 

Let us quantitatively discuss the range of applicability of the above formalism for real magnetic systems. To perform 
realistic calculations we use multiple scattering theory (see, for instance[8| and references therein) to obtain 'bare' J q . 
In this theory a key equation is 

r(s) = [P(e)-S}- 1 (14) 

where r(e) is scattering path operator, P(s) is the inverse of the one-site scattering matrix and S is the matrix of 
structure constants. It is convenient to separate spin structure explicitly: 

r(e) = T (e) + T(e)<r,P(e) = p(e) + p(e)«7. (15) 

Then the total energy variation with respect to the deviation of magnetic field at the site i can be presented in the 
static linear response scheme as[|| 

SE = ~J F delmTiL {Sp T m } , (16) 
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where Too is a matrix element of vector component of the full scattering matrix at site 0. Using the sum rule for r 
matrix introduced in Ref.Q, we can obtain the following relation for the 'bare' exchange coupling 

-l 

( Tt - Ti )oo= ( 17 ) 
i- J" cfalmTr, (TT - 1% [x], 1 (T* - T% , 

This expression properly takes into account the different energy dispersion for different spins in band magnets. Ea ll7l 
can be directly generalized for the non-collinear or spin spiral ordering. In this case general 4x4 supermatrix x should 
be build and inverted. Eg 1 171 can be seen as a generalization of long- wave approximation result [6j for the magnets 
with arbitrary degree of spin localization. Below we used the local density approximation with a linear muffin-tin 
orbital technique (LMTO) in the atomic sphere approximation ^4j to calculate the effective exchange according to 
EdlTl the SW spectrum 

= [Jq ~ ^0] /HI, (18) 

and the SW spectrum of the localized model as 

< = [4 - A] /m, (19) 

with long-wave J q determined in Ref.jl @. Eq|T7| is a static analog of EqQ] implemented in multiple scattering 
technique. Both the tetrahedron scheme and fast Fourier transform methods have been used for the Brillouin zone 
integration^. The matrix form of Eg 1 171 and a proper symmetrization of the exchange matrix from Eg. l|ll|) are 
essential ingredients of such calculations. The computation of the inverse of w q with only s,p and d basis set is 
numerically more stable compared to Xq\ but it is still affected by the type of approximation used for the effective 
tight-binding LMTO Hamiltonian (see below) . Three FM systems with entirely different degree of localization of the 
local moment were considered: Gd, Ni and Fe. An important observation is that in Gd (highly localized moments, 
small SW dispersion) the approximation based on the assumption A = Tr (x q — X ) Xo 1 << 1 i s completely fulfilled 
(A at the zone boundary is less then 0.01), whereas in Ni (moderately localized moments, large spin wave dispersion) 
it is not valid at all (A (Y) ~ 0.6), so that the corresponding matrix estimations using Eo. (|10|) demonstrate a strong 
enhancement of SW spectra at larger q (about 50% at Y point with two center approximation of LMTO 01 and 
70% using third order Hamiltonian with combined correction terms) . Fe is an intermediate case where the maximum 
of A is 0.30. This result indicates that the previous spin spiral calculations, where the derivative with respect to 
exchange-correlation field instead of local moment was used, the dispersion at finite q was underestimated. However, 
in the small q regime (SW stiffness), as we have shown above (see Eo llO|) . the old results are perfectly correct. 

It is important to analyze the effective exchange coupling between atoms in real space. Our results indicated that 
in Fe and Ni a main contribution coming from the renormalization of the first nearest neighbour exchange, so that in 
bec Fe it is increased from 7^ = 16.6 meV to Joi = 19.4 meV, whereas in fee Ni Jq 1 =2.7 meV to Joi=8.3 meV. 

Such a clear difference in the results indicate that the removal of the long-wave approximation can serve as an 
indicator of the degree of localization (parameter (x q — Xo) Xo~ la k° ve ) m different metallic magnets. In Fe and Ni, 
for instance, it at least partially explains why long wave mean field (MF) estimates predict such a small T c in FM Ni: 
300-350K in Ref.|4l IsL Il4j|. 340 K in present calculations, while experimental result is 630 K. First of all, long- wave 
approximation is suitable for such a 'localized' system as Fe, and the corresponding change in the nearest neigbor Joi 
is relatively small (correspondingly the increase of T c is expected to be small). FM Ni represents a rather itinerant 
system and any local approach (long wave approximation in particular) might produce a large error. In our case, a 
large increase in Joi for Ni indicates that the traditional MF approach (or any other approach which is based on 'no 
short-range order' assumption) is not applicable for the itinerant systems in general, predicting very high T c (above 
1000K). Such a number is also not consistent with the local density approach because latter does not allow to have 
T c in fee Ni larger then 500-540K (600-640K using gradient corrections) and this number can be considered only as 
an indicator of non-applicability of MF approach. 

So, the formalism above (see Eg. l|17[) ) tends to raise T c by inhibiting short-wavelength fluctuations and thus in- 
creasing the short-range order. Since this estimate already includes classically all wavelengths (except longitudinal) in 
the MF approximation and the final T c in Ni is expected to be very high, we have thus far first indication that other 
important mechanisms (with much 'softer' modes) are still missing in the modern finite-temperature magnetism theo- 
ries. Our preliminary SD simulations ^lj revealed the importance of the short range order effects in Ni and we believe 
that the proper definition of pairwise and non-pairwise interactions, which include both short and long range wave 



i- I'&ImTr, (2*-!% [j 



^T k T ^ +q 
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lengths correctly, is vitally important for the reliable description of the itinerant magnetism at finite temperatures in 
3d metals. 

In summary, we presented a technique for the calculation of the effective exchange parameters without any assump- 
tions about the form of magnetization density and the degree of magnetic order. Our results for elementary magnets 
(Fe, Ni and Gd) indicated that the most significant improvement is obtained for exchange coupling between nearest 
magnetic atoms and for SW spectrum at finite wave vectors. The large error of previous long- wave approximation 
is determined for fee Ni. This result indicates also that some important effects (which are essential for the correct 
description of T c in itinerant magnets) are still missing from the current magnetism theories when the nearest neigh- 
bors Jij parameters determined at T=0 K are being used at the finite temperatures. In Fe and especially in Gd the 
proposed technique gives the result similar to localized model result, as expected. 

Overall, we expect that the emphasized above approach for the exchange coupling and SW spectra Ea. (|8lllll7|) 
will significantly improve the description of the nearest neighbors interaction or large q part of SW spectra in any 
magnet with not fully localized magnetism (m < 2/i B and metallic r s ). We expect an increasing importance of this 
formalism at the finite temperatures and in systems with strong magnetic short range order. 

The author thank N. E. Zein and K. D. Belashchenko for many stimulating discussions and P. Bruno for the 
opportunity to present the results of this paper at the seminar in his group in Halle, in May 2002. This manuscript 
has been authored by Iowa State University of Science and Technology under Contract No. W-7405-ENG-82 with the 
U.S. Department of Energy. 
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